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分 数 模 永 磁 电机 的 多 物理 场 耦合 分 析 


杨 凯 颜 建 虎 张 尚 坤 
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摘要 : 由 于 电机 的 运行 状态 与 特性 与 其 电 、 磁 、 热 等 物理 场 直 接 相 关 ， 因 此 对 电 
机 进行 多 物理 场 耦合 分 析 至 关 重 要 。 本 文 利用 三 维 有 限 元 法 对 分 数 槽 永 磁 电机 在 运行 
时 的 磁场 分 布 进行 分 析 计算 ， 并 由 磁场 分 析 结 果 计 算出 各 部 分 损耗 情况 ， 将 各 部 分 损 
耗 作为 电机 运行 时 的 热源 ， 对 电机 的 温度 场 进行 仿真 分 析 ， 得 到 电机 的 温度 分 布 情况 。 
并 在 此 基础 上 ， 对 一 台 样 机 进行 实验 分 析 ， 实 验 结果 与 仿真 结果 一 致 。 从 而 验证 了 对 
该 电机 的 多 物理 场 耦合 分 析 方法 的 有 效 性 ， 为 后 续 的 电机 结构 优化 黄 定 了 理论 基础 。 
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Abstract: Due to the running status and characteristics of the electric machines 


are directly related to its electrical, magnetic, thermal and other physical fields, multiple 


physical field coupling analysis is of great importance. In this paper, 3D finite element 


method is used to research on magnetic field distribution when the fractional slot 
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permanent magnet machine running, and get loss values by magnetic field analysis. Taking 


控制 。 temperature field of permanent magnet machine and got the temperature distribution of the 
machine. Then a prototype machine is used to experiment based on this. The validity of 
the multiple physical field coupling analysis method is confirmed through the comparison 
of the experimental results with the simulation results, laid a theoretical basis for the 
subsequent electric machine structure optimization. 

Keywords: Permanent magnet machine, fractional slot, loss, multi-physics field, 


finite element method 


司 家 自然 科学 基金 项 目 (51407094) 和 江苏 省 自然 科学 基金 项 目 (BK20140785) . 


201903.00111v1 


chinaXiv 


1 引言 


与 传统 的 电 励 磁 电机 相 比 ， 稀 土 永 磁 电 机 具有 
结构 简单 、 运 行 可 靠 、 体 积 小 、 质 量 轻 、 损 耗 小 、 
效率 高 、 电 机 的 形状 和 尺寸 灵活 多 样 等 优点 ， 因 而 
广泛 应 用 于 航空 航天 、 国 防 、 工 农业 生产 和 日 常生 
活 等 领域 。 对 电机 磁场 的 分 析 和 计算 是 电机 设计 的 
重要 环节 ， 然 而 ， 永 磁 电 机 在 运行 过 程 中 的 温 升 对 
运行 状态 产生 较 大 影响 ， 使 永 磁 体 过 热 而 产生 不 可 
逆 退 磁 。 所 以 在 永 磁 电机 设计 和 分 析 中 对 其 热 场 的 
分 析 和 计算 具有 重要 的 意义 "1。 

目前 对 电机 物理 场 的 计算 和 分 析 主 要 采用 有 限 
元 法 ， 分 析 过 程 易 于 实现 标准 化 ， 可 得 到 通用 的 计 
算 程 序 且 有 较 高 的 计算 精度 ， 因 此 特别 适合 求解 电 
机 这 类 边界 形状 复杂 、 存 在 材料 非 线性 的 磁场 问题 。 
要 计算 电机 内 的 温度 场 ， 首 先 应 准确 地 计算 电机 内 
各 部 分 的 损耗 ， 确 定 电机 的 热源 。 永 磁 电机 损耗 主 
要 包括 铁 耗 、 定 子 绕组 钢 耗 、 机 械 损 耗 、 转 子 的 涡 
流 损 耗 以 及 未 计 入 上 述 损 耗 的 其 他 附加 损耗 ， 这 
些 损耗 会 转变 成 热能 并 经 电机 冷却 系统 带 走 ， 同 
时 使 电机 各 部 件 温度 升 高 。 因 此 对 永 磁 电 机 损耗 的 
精确 计算 对 准备 计算 和 分 析 电 机 温度 场 有 着 重要 意 
义 呈 1。 

目前 ， 国 内 外 采用 的 电机 热 计 算 方 法 主要 有 简 
化 公式 法 、 热 网 络 法 、 温 度 场 法 、 参 数 辨 识 法 和 基 
于 运行 状态 参量 的 电机 温度 模型 ~ ”"。 哈 尔 滨 理 工大 
学 的 李 伟 力 等 人 对 大 型 电机 定子 三 维 流体 场 及 其 对 
温度 场 的 影响 和 大 型 同步 发 电机 定子 多 元 流 场 与 表 
面 散热 系数 之 间 的 关系 进行 了 数值 计算 与 分 析 ， 并 
将 电机 内 的 流体 场 与 温度 场 进行 而 合计 算 处 。 沈 阳 
工业 大 学 的 王 凤翔 等 人 基于 磁场 有 限 元 和 3D 流体 
场 分 析 对 高 速水 磁 电 机 进行 了 损耗 计算 ， 并 基于 流 
固 炮 合 分 析 对 电机 的 温 升 进行 了 计算 中。 沈阳 工业 
大 学 的 唐 任 远 等 人 采用 三 维 有 限 元 法 对 永 磁 风力 发 
电机 进行 了 三 维 温度 场 的 计算 与 分 析 ， 验 证 了 采取 
热源 、 绕 组 等 效 导热 系数 的 计算 方法 的 可 行 性 '"。 
悉尼 科技 大 学 的 D. G. Dorrell 等 人 采用 2D 有 限 元 
分 析 方 法 研究 了 高 速 无 刷 永 磁 电机 线圈 中 的 趋 肤 效 
应 和 临近 效应 5。 瑞典 皇家 理工 学 院 的 0. Aglen 等 
人 对 高 速 发 电机 进行 了 损耗 计算 和 热 分 析 ， 结 果 表 
明 在 转子 中 加 入 铜 网 可 以 降低 损耗 值 "”。 芬 兰 阿尔 
托 大 学 的 Z. Kolondzovski 采用 2D-3D 数值 相 结 合 
的 方法 对 高 速 永 磁 电 机 进行 了 热 分 析 并 采用 有 限 元 

转子 动力 学 建 模 的 方法 对 电机 进行 了 机 械 性 能 的 分 
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析 "1。 

本 文 以 一 台 额 定 转速 为 1 500r/min 的 4 极 15 槽 
的 永 磁 同 步 电 机 为 例 ， 采 用 三 维 有 限 元 法 对 电机 的 
磁场 、 损 耗 和 温度 场 进 行 计算 ， 并 对 电机 在 稳定 运 
行 时 的 磁场 和 温度 场 进行 了 分 析 。 在 此 基础 上 与 样 
机 进行 了 实验 对 比 ， 为 后 续 电 机 的 结构 优化 设计 和 
分 析 黄 定 了 理论 基础 。 


2 永 磁 电机 的 磁场 分 析 


对 于 永 磁 电机 来 说 ， 其 磁场 的 准确 计算 分 析 是 
研究 电机 性 能 的 重要 环 市 ， 通 常 采 用 磁 路 法 对 电机 
进行 简化 计算 ， 然 而 计算 精度 较 差 ， 且 无 法 解决 电 
机 内 部 所 有 的 问题 。 因 此 ， 在 进行 永 磁 电 机 磁场 分 
析 时 ,通常 采 用 有 限 元 法 来 进行 分 析 和 和 计算。 考虑 
到 对 电机 损耗 的 准确 分 析 ， 三 维 有 限 元 法 成 为 主要 
的 计算 和 分 析 方 式 ， 用 三 维 有 限 元 分 析 磁 场 和 涡流 
场 ， 主 要 基于 4-9 法 ， 其 基本 方程 为 


WD nd (1) 
式 中 ，4 为 磁 矢 量 位 :9 为 电 标量 位 ;0 为 电导 率 ，; 
v 为 磁 阻 率 ， 为 线圈 的 电流 密度 。 
3 ” 永 磁 电 机 温度 场 分 析 


3.1 温度 场 理 论 分 析 
根据 传 热学 的 基本 理论 ， 在 直角 坐标 系 下 ， 电 
机 内 的 稳 态 温度 场 求解 可 归结 为 如 下 的 边 值 问 题 “" 


0 FR .0 站 a pe ee 
oOx\ x) 0 作 0 oz 0z 


-一 =0 (2) 


-4 -a(T -1) 
on 


式 中 , X44,, 4, 分 别 为 x,， y, z 方 向 上 的 导热 系 

数 ，gv 为 热源 密度 ， 是 电机 各 部 件 单位 体积 产生 的 

损耗 ，c 为 对 流散 热 系 数 ，7 为 周围 流体 的 温度 。 
由 变 分 原理 可 知 ， 式 (1) 的 等 价 变 分 方程 为 


ary for /ory 
加 + 加 必 | 台 | 六 
六 Ox Oy 0z 


CD- 了 


f Tav dy + or -27:7)dS =min 
V S3 


201903.00111v1 


chinaXiv 


对 其 作 离 散 化 处 理 后 ， 可 得 到 三 维 温度 场 有 限 
元 方程 


KT=F (4) 


式 中 , 了 为 求解 域内 全 部 市 点 温度 所 形成 的 温度 阵 
列 ; 到 和 五 分 别 为 总 体系 数 和 矩阵 和 总 体 右 端 矢 量 。 
对 方程 进行 求解 后 可 得 到 求解 域内 全 部 市 点 的 
温度 值 。 
3.2 电机 热源 的 分 析 
定子 铁心 的 热源 主要 来 自 于 绕组 铜 耗 和 铁心 损 
耗 "1。 通 过 电磁 计算 可 以 求 得 电机 铁心 中 的 磁 
场 分 布 ， 其 铁 耗 的 计算 公式 为 


nelem 


Wos = 之 Wi(B,f) + Wai(B, 7)) (5) 


式 中 ，W,(B) 为 单元 磁 琵 损耗 ， 丈 (4B) 单元 焦耳 损耗 。 

由 式 (5) 可 以 看 出 ， 在 分 析 铁 耗 时 ， 将 铁 耗 
分 为 磁 请 损耗 和 焦耳 损耗 两 部 分 ""。 铁 心 各 部 分 磁 
通 密度 的 准确 度 对 于 铁 耗 计算 的 准确 性 有 重要 影响 。 
由 于 永 磁体 上 的 涡流 损耗 相 较 于 电机 的 其 他 损耗 较 
小 ， 因 而 在 计算 中 忽略 ， 同 时 忽略 不 计 转 子 铁心 中 
产生 的 铁心 损耗 。 
3.3 电机 各 部 分 导热 系数 及 边界 条 件 

电机 中 的 热量 传递 方式 主要 有 热传导 和 对 流 换 
热 ， 即 导热 和 对 流 的 综合 过 程 ，""。 由 传 热 的 基础 
知识 可 知 ， 上 述 过 程 与 介质 的 导热 系数 和 表面 传 热 
系数 直接 有 关 。 由 传 里 叶 定 律 可 知 ， 导 热 系数 是 当 
温度 梯度 为 1 时， 单位 时 间 内 通过 单位 面积 的 导热 
量 。 导 热 系数 的 大 小 与 材料 的 性 质 有 关 ， 同 一 材料 
的 导热 系数 随 温 度 、 压 力 、 多 和 孔 性 和 均匀 性 等 因素 
而 变化 ， 其 中 温度 是 决定 性 因素 。 对 于 绝 大 多 数 物 
质 而 言 ， 材 料 温度 尚未 达到 融化 或 气 化 以 前 ， 导 热 
系数 可 以 近似 地 认为 是 线性 规律 变化 ， 即 : X46=Ao(l1+ 
20。 甚 中 ，?o 指 温度 为 零 时 的 导热 系数 ， 2 是 由 试 
验 确定 的 常数 。 本 文中 永 磁 电 机 主要 材料 的 导热 系 
数 见 表 1。 

计算 定 转 子 全 域 温 度 场 时 ， 需 要 对 气 际 进行 特 
殊 处 理 ， 可 以 用 热 方 式 换 热 的 效果 等 价 气 隙 中 对 流 
方式 换 热 的 效果 。 有 效 导热 系数 可 按 下 述 方法 求 取 。 

假定 定子 内 表面 和 转子 外 表面 为 光滑 圆柱 面 ， 
而 气 隙 中 的 雷诺 数 Re 可 表示 为 


Wyo 
Ne (6) 
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表 1 电机 主要 材料 导热 系数 


Tab.l Thermal conductivity of main materials 


部 件 名 称 导热 系数 /(W/m : K) 
定子 铁心 39 
空气 0.0305 
定子 绕组 385 
水 磁体 9 
转子 极 靳 4.43 
转子 斩 部 45 


式 中 ，ws%i 为 转子 的 圆周 速度 ; 6 为 气 际 的 长 度 ; E 
为 空气 的 运动 粘度 。 
临界 雷诺 数 Re 的 表达 式 为 


Re, 42 人 (7) 


式 中 ，R; 为 定子 内 径 。 
当 Re < Re 时 ， 气 隙 中 的 空气 流动 为 层 流 ， 
有 效 导 热 系数 her 等 于 空气 导热 系数 ; 当 Re > Re.. 
时 ， 气 隐 中 的 空气 流动 为 率 流 。 计 算得 出 
Re =1325 > Re =839.2 (8) 


此 时 ， 气 际 中 的 空气 流动 为 大 流 ， 气 隙 的 有 效 
导热 系数 为 


A 生 汉 0.00197 > Re m3336) (9 ) 


ee 


式 中 ,7 = rwyRi; 16 为 转子 外 径 ，Ri 为 定子 内 径 。 
通过 计算 得 出 气 险 有 效 导 热 系 数 为 0.103(W/ 
mK)。 此 外 ， 电 机 温度 场 计 算 的 边界 条 件 包括 : 
生 热 边界 、 散 热 边界 和 绝缘 边界 。 


4 ”模型 建立 及 网 格 剖 分 


电磁 场 分 析 的 基本 流程 ， 可 分 为 前 处 理 、 解 
算 和 后 处 理 三 个 阶段 。 前 处 理 阶 段 主要 包括 建立 分 
析 模 型 (电机 磁场 有 限 元 模型 和 功率 变换 器 电路 模 
型 )、 设 定 材料 属性 、 边 界 条 件 和 仿真 参数 以 及 网 格 
剖 分 等 儿 个 步骤 。 在 后 处 理 阶段 ， 可 以 将 电机 磁场 
有 限 元 模型 中 每 一 个 单元 或 节点 的 磁场 数据 以 及 电 
路 模型 中 的 电流 、 电 压 等 数据 提取 出 来 进行 数据 分 
析 和 处 理 。 

本 文中 所 研究 的 分 数 模 永 磁 电 机 模型 和 网 格 
剂 分 分 别 如 图 1 和 图 2 所 示 。 图 3 为 电机 的 外 电路 


模型 。 
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图 1 电机 模型 
Fig.l Model of machine 


图 2 网 格 剖 分 
Fig.2 Mesh of machine 


FEM FEM FEM 
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图 3 外 电路 模型 
Fig.3 Circuit model 


5 电机 磁场 仿真 结果 


电机 空 载 电动 势 及 齿 槽 转 和 矩 波 形 分 别 如 图 4 和 
图 5 所 示 。 由 于 采用 开口 醒 设 计 ， 其 空 载 电动 势 具 
有 一 定 程度 的 耸 谐 波 。 齿 槽 转 矩 是 永 磁体 不 通电 时 
永 磁体 和 铁心 之 间 相 互 作用 产生 的 转 矩 ， 是 由 永 磁 
体 与 电 枢 齿 相 互 作用 力 的 切 向 分 量 的 波动 引起 的 。 
齿 槽 转 矩 导致 转 和 矩 波 动 ， 影 响 系统 的 控制 精度 ， 同 
时 产生 振动 和 噪声 站。 
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一 一 0 相 一 一 Vv 相 RS wW 相 


空 载 电动 势 /V 


0 30 60 90 120 150 180 
机 械 角度 /(*) 


图 4 ” 空 载 电动 势 


Fig.4 No-load electromotive force 


齿 槽 转 矩 /(N*m) 
= 
尼 一 ID 


-0.1 
-02 
-0.3 
0 30 6 9% 120 150 180 
机 械 角 度 /() 
图 5 齿 权 转 矩 


Fig.3 Cogging torque 


图 6 为 电机 磁 密 的 整体 分 布 情况 ， 可 以 看 出 转 
子 上 有 部 分 饱和 现象 ， 其 中 定子 铁心 的 平均 磁 密 为 
0.87T， 转 子 铁心 平均 磁 密 为 0.95T。 图 7 为 电机 定 
转子 磁 矢 量 图 ， 通 过 图 7 可 以 分 析 电 机 漏 磁 分 布 情 
况 。 通 过 以 上 对 电机 磁场 的 计算 分 析 可 以 得 出 ， 电 
机 的 铜 耗 和 铁 耗 分 别 为 34W 和 5.77W。 


图 6 磁 通 密度 云图 
Fig.6 Magnetic flux density 


图 8 为 电机 额定 负载 时 的 输出 转 矩 。 此 时 ， 
电机 的 输出 转 矩 为 3.3N .mm， 但 是 存在 转 矩 脉动 ， 
该 脉动 是 由 于 齿 模 转 矩 和 控制 方式 产生 的 ， 约 为 
7%, 


201903.00111v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


7 ” 磁 通 密度 矢量 图 


Fig.7 Flux density vector 


负载 转 矩 (Nm) 


0 30 60 90 120 150 180 


机 械 角 度 /(”) 
8 负载 转 矩 
Fig.8 Load torque 


6 永 磁 电 机 温度 场 计算 结果 分 析 
通过 磁场 分 析 得 到 的 电机 损耗 作为 热源 对 其 温 


度 场 进行 仿真 计算 ,结果 如 图 9 和 


图 10 所 示 。 从 仿 


真 结果 中 可 以 看 出 ， 电 机 的 最 高 温度 出 现在 定子 此 
部 ， 大 约 为 83.8C 。 这 是 因为 绕组 绝缘 及 绝缘 介质 
的 导热 系数 很 小 ， 而 且 线圈 绕组 也 是 一 个 大 功率 的 
发 热源 ”"。 有 具体 各 部 分 温 升 情况 在 表 2 中 给 出 。 


9 ”定子 铁心 及 绕组 温度 分 布 


Fig.9 The temperature distribution of stator core 


and winding 


图 10 转子 及 永 磁体 温度 分 布 


Fig.10 The temperature distribution of rotor and 


permanent magnet 


表 2 电机 温度 场 计 算 值 


Tab.2 Thermal field calculation results 


区 域 温 升 /C 

最 高 温 62.47 
定子 绕组 最 高 温 升 62.47 
定子 齿 部 平均 温 升 62.34 
定子 斩 部 平均 温 升 61.85 
永 磁体 最 高 温 升 41.66 
永 磁体 平均 温 升 41.25 
转子 罗 部 平均 温 升 40.60 


电机 运行 时 发 热 均 来 自 甚 损耗。 其 中 定子 和 绕 
组 既是 发 热 部 件 也 是 传 热 部 件 ， 别 的 部 件 为 传 热 音 
件 。 随 着 电机 温度 的 升 高 ， 它 与 周围 介质 之 间 的 温 
差 逐 渐 增 大 ， 散 失 到 周围 空间 中 的 热量 相应 地 也 逐 
渐 增多 ， 本 身 温度 升 高 的 速度 则 逐渐 变 慢 ， 最 后 电 
机 产生 的 全 部 热量 都 传 给 周围 介质 ， 达 到 热 稳定 状 
态 。 通 常 每 小 时 温度 变化 小 于 1C 时 ， 即 认为 已 经 
达到 热 稳定 状态 。 从 图 11 的 电机 定子 温度 变化 仿真 
曲线 可 以 看 出 ， 电 机 运行 一 段 时 间 后 ,温度 最 大 值 
达到 了 847C 。 


100 
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时 间 /s 
11 电机 定子 温度 变化 仿真 曲线 和 实验 数据 


Fig.ll Simulation curve of motor temperature and 


温度 /C 


experimental data 
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图 12 给 出 了 电机 样机 模型 ， 通 过 对 样机 进行 实 
验 得 出 的 定子 铁心 温度 测量 数据 可 知 ， 当 电机 运行 
70min 后 ， 定 子 铁心 中 的 平均 温度 达到 了 88C ， 之 
后 温度 基本 保持 恒定 ， 与 计算 结果 基本 一 致 。 


7 a 


(a) 定子 模型 
图 12 样机 模型 
Fig.12 Prototype motor 


(b) 转子 模型 


本 文通 过 三 维 有 限 元 法 对 分 数 槽 永 磁 电 机 的 磁 
场 进行 了 分 析 与 计算 ， 并 基于 磁场 计算 推导 得 出 电 
机 损耗 分 布 ， 利 用 各 部 分 损耗 作为 热源 对 电机 的 磁 
热 耦 合 展开 分 析 。 在 此 基础 上 ， 对 样机 进行 实验 测 
量 ， 从 实验 结果 看 出 ， 电 机 模型 各 部 分 的 温度 分 布 
趋势 与 计算 数据 基本 一 致 ， 为 后 续 电机 的 结构 优化 
设计 黄 定 了 理论 基础 。 
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